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. Abstract 

m 

In the stochastic limit the resonances play a fundamental role because 
they determine the generalized susceptivities which are the building blocks 
\ of all the physical information which survives in this limit. There are two 

■ sources of possible divergences: one related to the singularities of the form 

factor; another to the chaoticity of the spectrum. The situation will be 
, illustrated starting from the example of the discrete part of the hydrogen 

atom in interaction with the electromagnetic field. 

^ ■ 1 Introduction 

^ , . Beyond the scheme microscopic-mesoscopic-macroscopic, there are many levels 

^ '' of description, probably an infinite hierarchy, in which the specific properties of a 

given level express some kind of cumulative or collective behaviour of properties 
of systems corresponding to lower levels. 
. These cumulative phenomena are, typically, nonlinear effects. In absence 

' of generally applicable methods, one introduces asymptotic methods which ap- 

proximate the values of individual quantities of physical interest. Among these 
asymptotic methods, the scattering theory is concerned with the long-term be- 
haviour of physical system, for which t cx), and the perturbation theory is 
concerned with weak effects, for which A — > 0, where A is a parameter upon 
which the interaction Hamiltonian depends. 

The stochastic limit 1^ puts together the scattering and perturbation the- 
ories by studying the long-term cumulative effects of weak actions, working as 
a magnifying glass of all phenomena pertaining to the scales of magnitudes we 
are interested in and as a filter of those pertaining to all the remaining scales. 

*In Memory of I. Prigogine. Presented at the CFIF workshop on time asymmetric quantum 
theory: the theory of resonances, Lisbon, Portugal, 23-26 July, 2003. 
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One starts from the usual quantum Hamiltonian dynamics in interaction rep- 
resentation, see Sect|21 Depending such dynamics on a parameter A, one rescales 
some parameters in the associated evolution (typically time) and obtains in the 
limit A — > a new dynamical system driven by a singular Hamiltonian. The 
new unitary evolution is an approximation of the original one which preserves 
much nontrivial information on the original complex system, see Sect|31 

The limit Hamiltonian is a functional of some white noise, see SectJS] The 
idea is that, if we look at the fast degrees of freedom of a nonlinear system with 
a clock, adapted to the slow ones, then the former look like an independent 
increment process, typically white noise. 

Moreover, in the limit we shall have not a single, but an infinity of indepen- 
dent quantum noises, one for each Bohr frequency of the system. This is the 
stochastic resonance principle, see Theorem 0] 

A standard scheme to describe dissipation and irreversibility passes through 
the so-called open system (or system-reservoir) approach, the basic physical idea 
of which is that dissipation an irreversible behaviour arises when two systems, 
traditionally denoted S {system: slow degrees of freedom) and R {reservoir, fast 
degrees), interact and one of them (i?) exhibits some macroscopic or chaotic 
features, see Sect0] 

In the open-system approach the master equation corresponds to the adia- 
batic elimination of the fastly relaxing variables, a technique also called coarse 
graining, obtained just by taking the partial expectation of Langevin equations 
(stochastic limit of the Heiscnberg evolution) with respect to the reference state 
of the master field. 

The stochastic limit goes far beyond the master equation because it does not 
eliminate the fast degrees of freedom. This allows to estimate the probabilities of 
some collective states, or more generally, the behaviour of a complex (nonlinear) 
system with many degrees of freedom in terms of relatively few functions of the 
microscopic characteristics of the quickly relaxing degrees of freedom. According 
to the interpretation, these functions are called order parameters, kineticies, 
susceptibilities or susceptivities, transport coefficients, etc, see Sect|31 

There is one generalized susceptivity factor for each Bohr frequency; its real 
part is a (5-function and its imaginary part a generalized Hilbert transform, both 
over each resonant surface, see Eq. (|31|l . We shall establish some connections 
between these real and imaginary parts in Sect El 

The transport coefficients appear as Ito correction terms in the white noise 
Hamiltonian equations, which introduce quantum mechanical fluctuation-dissi- 
pation phenomena in the evolution. The imaginary part of the Ito correction 
term corresponds to a global shift in the spectrum of the system Hamiltonian, 
whereas the evolution corresponding to the real part of the Ito term is a con- 
traction, i.e. dissipative. 

Finally, in Sect [7| we shall illustrate the situation with a concrete physical 
example. We will study the bound states of the hydrogen atom in interaction 
with the electromagnetic field. In this case we shall obtain some sufficient 
conditions for the existence for the generalized susceptivity factors and calculate 
them explicitly. 
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2 The Interaction Representation 

Let H he a, Hamiltonian and let ^p(t) be a solution of the corresponding Schrodinger 
equation (?i = 1) 

Given a decomposition of H into a free and an interacting parts, 

H = Ho + Hi, 

the wave function ^/(i) in the interaction representation is defined by 

^j{t) e'*""iP{t) = e''""e~""-"°+"'^o, (1) 

and it satisfies the Schrodinger eq. in the interaction representation 

dtMt) = ~iHi{t)i,i{t), Hj{t) e^*^»i/,e-'*^». (2) 

The solution of Eq. (O with initial condition V'/(^o) is given by tpi{t) = 
U{t,to)ipi{tQ), where the propagator U is 

U{t, to) = e**-f^og-»(t-to)(Ho+if7)e"rfoffo^ (3) 

which verifies 

Uit,to) = U{t,s)U{s,to), U{tM* ^U[to,t), U{t,t)^I, (4) 

U{t + r, to + r) = e^'^^°C/(tito)e-^'^^'' (5) 

dtU{tM) = -iHi{t)U{t,ta), U{to,to) = I. (6) 
The S matrix is formally defined as the limit 

S= lim U{t,-t). (7) 

In interaction representation the Heisenberg evolution is 
XtoiA)^At ■.= U{t,to)*AU{t,to) 
and satisfies the flow equation 

jutoiA)=js,toJtA^), (8) 

whose differential form is the Heisenberg equation in the interaction representa- 
tion 

dtAt = -t[Hiit),At]. (9) 

In what follows, by a dynamical system we mean a pair {7i, H}, where Ti. is 
a Hilbert space and H a Hamiltonian, or we mean a triple system 

{H,Ho,U{t,s)}, 

where {Ti., Hq} is an integrable dynamical system and U{t, s) is a propagator 
satisfying the cocycle equations and If U{t,s) is differentiable with 

respect to t, then the two definitions are equivalent. However, the generalized 
definition also allows the possibility for U(t, s) to satisfy a stochastic differential 
or white noise Hamiltonian equation. 
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3 A First Approach to the Stochastic Limit 



The stochastic hmit approximates the fundamental laws themselves in the fol- 
lowing sense: One starts with a family of quantum dynamical systems {Hx , Hq^^ , J//'^'' }, 
depending on a parameter A, verifying the symbolic relation^ 

liin{nx,H^^\ul^l} = {H, Ho, C/J, (10) 

where {Ti., Hq, J7t} is a new quantum dynamical system. 

The equation satisfied by the dynamical evolution Jj/^-* is^ 

dtU^^^ = -tH^,^\t)U^^\ Ul^^^^l, (11) 

interpreted as a Schrodinger equation in interaction picture (and not as the usual 
Schrodinger equation with time-dependent Hamiltonian) . With the change of 
variables t^X^t Eq. takes the form 

dtU[% = -^y^Hf^\t/\^)ul%, C/(^) = /. (12) 

The first step of the stochastic limit is to prove that 

hm^i/f)(VA2) = i7,, (13) 

this implies u'^'p^2 Ut and Ut is the solution of 

dtUt = -iHtUt. 

The following lemma says us that when Hj'^\t) is independent of A and an 
integrable function of t, then the limit Hamiltonian Ht is simply a multiple of 
the (5-function in t characterizing the white noise correlations. 

Lemma 1 For any integrable function F onM. one has 

limi,F(l^)=^(r-t)/^F(.)da, (14) 

in the sense that, for any bounded continuous function -0 on M, 

lun ^ ^F (^^^ V^(r) dT = ^(r) ^ da. (15) 

Why the rescaling t — > t/X^ as the new time scale is best understood by 
considering second-order perturbation theory: 

^In Sect|5]wc will sec that the limit 1101 has to be understood in the sense of correlators. 
■^Note that uj:'^^ is the adjoint of the backward wave operator at time t, i.e. f/^^' =; 
Q5^^*(t) — > fi*^'*, where C^^^ is the backward wave operator. 
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Lemma 2 Denote by (•) = ($,•$) the expectation value with respect to a fixed 
vector Suppose that the Hamiltonian has mean zero, is time-translation in- 
variant and such that the function s i— > {Hj{0)Hj{s)) is integrable, i.e. 

{Hjit)) = 0, 

{Hl{ti + s) ■ --Hiitn + s)) = (Hl{ti) ■ ■■Hl{tn)), 

[ \{HiiO)Hi{t))\dt<cx,. 

JR 

Then the expectation value of the second- order term of the iterated series for 
U^'\ ^.e. 

/" dt^ r dt2{Hi{t^)Hi{t2)), (16) 
JVi Jo 

has a finite nonzero limit as A — > and t ^ oo iff 

lim A^i = r = c*VO(<oo). 

In this case the limit of fJ6|) as A — > is equal to 

-T f ds{Hi{0)Hi{s)). (17) 

The integral in Eq. (|17|l can be interpreted as a quantum transport coeffi- 
cient. In several cases (e.g. in the weak coupling, but not in the low-density) 
the transport coefficient allows us to compute the lifetimes and energy shifts in 
agreement with the Fermi golden rule. 



4 Open Systems 

Given two quantum dynamical systems, the system S — {Ti.s,Hs} and the 
reservoir R =^ {Hr, Hjj}, the quantum dynamical composite system will be of 
the form 

{Hs «) Hr, Hsr ^Hs(E>1r + 1s(^Hr + Hi}, 

where the interaction Hamiltonian Hi contains all the new physics, with respect 
to the isolated systems, while Hq — Hs + Hr (resp. Hs, Hr) is the free 
Hamiltonian (resp. of 5*, R). 

As reservoir R we will consider a quantum field Ok, a^t over M'', d > 3, for 
which the quantities of physical interest are the correlators (a^'^ • • • a^" ) , where 
(•) is a expectation value or state (i.e. a positive linear functional on the field 
algebra). Giving an expectation value is equivalent to giving a representation 
of the field algebra in a Hilbert space H and a unit vector $ G 7i such that 
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Recall that the field Ofc, a^, together with the expectation value (•) is called 
a mean zero Gaussian field if 

{a'j}_^ ■ ■ ■ a^" ) = 0, if n is odd, 
K\ • • • ) = J2 , ii ; . . . V, Jp) {a;i^ al7^ ) • • • {a^^^ a^J^ ) , 

•Po(2p) 

where 7^o(2p) is a subset of the ordered partitions (ii, ji; . . .ip, jp) of the set 
(l,...,2p) such that ii < ■■■ < ip and < ja, and e{ii, ji] . . .ip, jp) is a 
complex number. 

For a mean zero Gaussian field the state (•) is determined by its covariance 

The a\ are operator-valued distributions. The translation from the distribu- 
tion to the operator language requires integration against suitable test functions 
9- 

/ gik)a+dk, A{g)= / g{k)akdk. 
Jm Jm 

The rigorous meaning of the multiplication of distributions is the multipli- 
cation of the corresponding operators A'^{g). 

In what follows we shall assume that the reservoir i? is a mean zero Gaussian 
quantum field Ok, a'^ over M.'^, o? > 3, which satisfies the q-commutation relations 

[ofc, a^,]g = afcfl^, - qa'^,ak = 6{k - k') 

for some complex number q. In this case are called annihilators and a^, 
creators. We speak of a boson field if g = 1, of a Fermi field if g = —1, of a 
Boltzmann field if g = 0. 

If ak,a^ are boson creation and annihilation operators, there exists a real- 
valued function ^{k) such that the associated free Hamiltonian has the form 

Ho = J uj{k)a+{k)a{k)dk, (18) 

in the sense that the commutator with Hq of any polynomial in the field opera- 
tors coincides with the commutator of the same polynomial with the right hand 
side of Eq. (|18|l . In particular, for all k we have 

dtakit) - -i[Ho, ak{t)] = -icokakit) ak{t) = £-"'"'=0^(0). 

The function uj{k) is called the free 1-particle Hamiltonian. Usually one wants 
uj{k) to be positive and vanishes nowhere. For example, fc^/(2m), \/k^ + rri^ , 
\k\. 

Consider a system-reservoir Hamiltonian 

HsR^ = Hs®Ir + Is®Hr + Hi, 
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for which the interaction Hamiltonian Hi is of dipole-type, i.e. it has the form 
Hj = [ dk{D{k)(E)a+{k)+D+{k)(»a{k)}, (19) 



where {D{k) : k G M } is a family of system operators (acting on TLs) called 
the response terms or currents and containing local information about the in- 
teraction. 

If the free system Hamiltonian Hs has a discrete spectrum, i.e. 

Hs ^^EnPn = ^£n\£n){£n\, 
n n 

the time evolved interaction Hamiltonian becomes 

Hi{t) = I dkY^ P,nD{k)Pn ® e**("'=+""-''")a+ + h.c. (20) 

rn,n 

Introducing the operators 

DUk):= PmDik)Pn= (en|i^(fc)|£m), |e™)(en|, (21) 

for which we have e**^s£)^g-»tffs = e'"'*'^ D^{k), Eq. ^ becomes 

Hi{t) ^ I dkY Du{k) ® e''^'^''-'^^a+ + h.c. (22) 

Moreover, assume that the system operators D^{k) verify the generalized 
dipole approximation^ i.e. 

D^{k) = Dt:^°'-{k) = g{k)D^, 

where 5 is a test function, the cutoff ov form factor describing the strength of the 
interaction of the system with the environment, and D^j is a fixed system oper- 
ator. Then, if there is only one Bohr frequency u>, the interaction Hamiltonian 
ifT^ becomes 



Hi^D+ ® A{g) +D®A+{g)^ dk {D+ ® g{k)ak +D® g{k)a+}. 
and its free evolution is 

Hi{t) = I dk {D+^gf,e-''^^^-^^ak+D^g,e'*^^--^^a+} ^^^^ 
D+at + Da+, 
where we omit the symbol (g) and 

at := J dfcgfce-**("'=-")afe. 
Let us show how quantum white noises arise as stochastic limits of free fields. 
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5 The Stochastic Resonance Principle 

A quantum field b^{t, k) G R, fc G K*^) with expectation value (•) is called a 
white noise if it is a mean zero Gaussian field with covariance 

{b'{t, ky'{t', k')) = 5{t ~ t')G,^,>{k, k'), 

where e'{k, k') is a positive definite distribution (as a function of the variables 

The convergence will be taken in the following sense. Let aj(fc) be a family 
of fields with cyclic vector parametrized by a real number A, and let a J be 
another field with cyclic vector $. a^(fc) is said to converge to in the sense 
of correlators^ as A — s- Aq if, for any natural n, for any fci . . . A:„ e R'', and any 
choice of ei . . . e„ e {0, 1}, one has, in the sense of distributions, 

lim ($A,alHfci)---4"(fcn)'fA) = •••4;$). 

A— »Ao 

It is well known that if, as A ^ 0, the correlators {^\,a'^{ki) ■ ■ ■ a^!^''{kn)^\) 
converge to a distribution F{ki . . . kn), then there exists a field {aj, 7i, $} such 
that F{ki ...kn) ^ ($, al\ ■■■all^). 

Coming back to the time evolved interaction Hamiltonian given in Eq. (|23|) . 
The Schrodinger equation in interaction representation takes the form, see Eq. 

dtU^^^ = -iH\^\t)U^^^ = ~iX{D+at + Da+)U[^\ (24) 

The rescaled evolution operator satisfies the reseated Schrodinger equation, see 
Eq. ini): 



(25) 



The next result shows that the rescaled fields jOt/A^ converge in the sense 
of correlators to a quantum white noise bf . The new fields operators b act on 
some new Hilbert space and are called the master fields or quantum noises. 

Theorem 3 If the field (k) is mean zero Gaussian and satisfies the q-deformed 
relations 

a{k)a+{k') - qa+{k')a{k) = 5{k - fc'), 
then the rescaled field 

\a^{t/\\k) :=ie±'"(^)*/^'a±(fc) 
A A 

converges in the sense of distribution correlators to a q-deformed white noise 
bf{k), i.e., satisfying 

btik)b+{k') - qb+{k')bt{k) = 27r 5(t - t') S{uj{k)) S{k - k'). 



•^The correlation functions are also called Wightman functions in quantum field theory or 
mixed moments in probability theory. 
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Thus, from Eq. 



1 



lim ^i/f ^ (t/X^) ^Ht = D+bt + Dh+ (26) 
and the white noise Hamiltonian equation is 

dtUt = -iHtUt = -i{D+ht + Dht)Ut. (27) 

One has to bring equation (|27|) to its normally ordered form, i.e. instead of 
the term btUt we would hke to have the term Ufbt- If the reference state of the 
a field is the Fock vacuum, then the normally ordered form of the white noise 
Hamiltonian (|27|l is the following equation 

dtUt = -i{D+Utbt + Db+Ut) - YUt, (28) 

where 

Y := j-D+D, (29) 

7_ := /" dt f dfce-^*("'=-'^)|g(fc)|2. (30) 

The operator Y defined by Eq. (|29|l . the coefficient of Ut in the term not 
including noise operators in Eq. (|28|l . is called the operator transport coeffi- 
cient, drift coefficient or Ito correction term. This term corresponds to a non- 
selfadjoint correction to the system Hamiltonian that can be considered as the 
prototype of the quantum mechanical fluctuation- dissipation relation: 

iHs iHs - -i-D+D = i{Hs + Im{j^)D+D) - Re{-f-)D+D. 

Since Hs commutes separately with the real and imaginary part, the imagi- 
nary part of the Ito correction term is a global shift in the spectrum of the system 
Hamiltonian. On the other hand, the evolution generated by —Re{'y^)D~^D is 
a contraction, in general nonunitary, i.e. dissipative. In fact, 2i?e(7_) = {btbf) 
is a measure of the strength of the fluctuations. 

For the generalized susceptivity factor 7_ defined by Eq. ()30|1 . using the 
identity 

' dt = = n5{uj) - iP.P.-, 



we can write 



ui — iO to 



TT{g, S{uJk - io)g) - i{g, P.P. ^ g) 

\9{k)? 



11 / dk\g{k)\^5{oJk-oj)-iP.'P. I 



(31) 



dk- 



Assuming that cj(fc) > almost everywhere, the argument of the (5- function 
in Re{'-f-) can be zero only if > 0. Since 2i?e(7_) = (btbf), the fact that the 
contribution relative to a given frequency vanishes means there is no master field 
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with that characteristic frequency. So, in the stochastic hmit only the master 
fields corresponding to positive Bohr frequencies survive. But Eq. ^ also 
shows that this is not the case for the imaginary part. In the stochastic limit 
the negative Bohr frequencies contribute with an overall red shift to the energy. 

For more than one Bohr frequency the time evolved interaction Hamiltonian 
is of the form, see Eq. (|^ : 

Hi{t) = f dkJ2 D^g{k)e''^'^''-'^^a+ + h.c. (32) 

This decomposition of the interaction Hamiltonian suggests that, before the 
limit, the original field splits into a family of effective fields, each of which is 
interacting with its own Bohr frequency lo and with the other effective fields. 

The following theorem says that, if it is possible to interchange the sum, the 
integral and the limit, the mutual interaction of these effective fields becomes 
negligible in the stochastic limit, and after the limit this mutual independence 
becomes exact. 

Theorem 4 Let a(t,k), a'^{t,k) be a mean zero Gaussian quantum field with 
respect to a given state (•). suppose that the process a'^{t,k) is stationary and 
its covariance matrix is time-integrable in the sense of distributions, i.e. 

<cx), f^geSiR"). 

For each real number lo define the new process 

a^{t,k) := e~**"'a(t,fc). 
Then the limit, in the sense of distribution correlators, 

lim Y flw {t/ , fc) = b^ {t,k), 

exists and is the white noise b^{t,k) of the same Gaussian type as a(t,k) and 
with covariance 

{hl{t,k)bi{t',k'))^5{t^t') [ dse'''''^{a^iO,k)a''is,k')). 

Moreover, the quantum white noises {b^^(t,k) : lo G M} are mutually indepen- 
dent. 

This illustrates the stochastic resonance principle: In the limit we shall have 
not a single, but an infinity of independent quantum noises, one for each Bohr 
frequency of the system. 

Thus, for each Bohr frequency we will have a generalized susceptivity factor 
given by Eq. H31I) . whose real part is a (5- function and its imaginary part is a 
generalized Hilbert transform, both over each resonant surface. Let us pass to 
study these real and imaginary parts. 



dt 



dkdk' f{k)g{k'){a%0, k)a'''{t, k')) 
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6 The Distribution {uj{k) — u)~ 



Let 



k e 



}{k) - Lo e 



be a C°°-function except perhaps for some closed set of singular points with 
d-dimensional Lebesgue measure zero. In this section we study, for each test 
function 4> G X'(M'^), the integral 



(jj{k) — uj 



.m) = 



m 



u}{k) — UJ 



and its Cauchy principal value or principal part 



(P.P. 



1 



uj{k) 



— ,(j){k)) = lim 



dk 



m 



\uj{k) — Lo\ >e 



uj{k) 



■ dk , 



just the imaginary part of the generalized susceptivity factor 7_, see Eg. (13111 . 
To this end, assume that Vtj(fc) ^ for each regular point k e R"^. Then 



the equation Lo{k) 



UJ 



defines a regular surface Sc of dimension d — 1 in 



a neighbourhood sufficiently small of each regular point verifying the equation 
and, in such neighbourhood, we can consider a local change of variables tp : 
(ui, . . . , Ud) (fci, . . . , kd) such that the surface Sc is given by fixing one of the 
tt's to the value of the constant c. So, for example, we can choose uj(k) — uj = ui 
and arbitrarily the U2, ■ ■ ■ ,Ud, but with the condition that the Jacobian J'ip{u) 
is different from zero. Then, by the change of variables theorem, we will have^ 



m 



■ dk 



Uj{k)-UJ J^-l{Rd} Ul 

being the last integral, by Fubini's theorem, equal to 



Jijjiu) du , 



1 

a Ul 



(t){ip{u)) J^{u) du2 ■ ■ ■ dud 



dui , 



(33) 



where \j}-'^{W^) = {ui x iIj-^ \ ui e [a, b]}. 

We can write the integral (|33|l in terms of differential forms. For it, consider 
the form fie of order d — 1 associated to the function W{k) = uj{k) — on the 
surface Sc by the equation 

dW KVlc^dv, (34) 

where dv = dki A • • • A dkd is the volume element in K*^. Such form fij, ex- 
ists in a certain d-dimensional domain containing the surface Sc because in a 
neighbourhood of any point of that surface one can introduce a local system 
of coordinates ui, . . . ,Ud such that one of these coordinates, for example Uj, is 

''We can assume that the support of {j> is contained in the neighbourhood where the local 
change of variables is given, else we consider a suitable locally finite partition of the unity. 
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the quantity W{k); then, if ip : (ui, . . . , Ud) ~* {ki, . . . , kd) is the corresponding 
change of coordinates, we wiU have 

dv — jTp{u) dui A • • • A duj^i A dW A rfuj+i A • • • A dud 

and then 

fie = A • • • A duj^i A cJmj+i A • • • A dud- (35) 

In particular, if in a neighbourhood of a given point we have dW{k)/dkj ^ 0, 
we can take as coordinates u 

ui ^ ki, . . . ,Uj ^ W, . . . ,Ud = kd , 

we have then 



Jip-^k) dW{k)/dkj 
and the form Vic defined in Eq. (|35|l becomes 

„ _ , ^ j_i dki A • • • A A dfcj+i A • ■ 



A dfcrf 



The form 0,^ verifying Eq. H34I) is not unique since we can add to it any 
form A orthogonal to dW , that is, such that dW A A = 0. Such forms A can be 
written as A = 7 A dW ^ where 7 is a certain form of order d — 2. 

We note also that the form fl^ does not depend on the choice of coordinates 
ui, . . . duj^i, Uj+i, . . . , Ud, but it does on the function W defining the surface Sc- 
For example, if instead of W(/c) = c one considers the equation Q!(fc) [Ty(A:)—c] = 
0, where a is a nowhere zero function, then d{a{yV — c)) = a dW + {W — c) da 
and along the surface Sc we will have 

= '^^ = 

" d{a{W-c)) a ' 

Coming back to the situation of formula H33|l . for coordinates ui = W, U2, . . . , Ud, 
by formula H35I) we have 

ilc — — J''Pi''^)\ui-c'^^^ A • • • A dud 
and from Eq. H33|) we obtain 



i{k) -uj ui 



dui . (36) 



In terms of the distribution S{W — ui) ~ S{uj{k) — lo ~ ui) formula (|36f) can 
also be written as^ 



" 1 



— {d(W-ui),(l)) dui 

Ui 



'For a definition of the distribution 5{W — ui) see, for example, the chapter 3 of |2]- 
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We note that the function 



$(ui)= / (/.(^(u))0„,, (37) 

that appears in formula is a test function belonging to ^'(K). Indeed, the 
change of variables ^p, being bicontinuous, transforms compact sets into compact 
sets, and vice versa. Moreover, since '0 is a C°°-function, (f) o ip & ViW^) and 
then its product with the C°°-function Jif) is also in ViW'-). Thus, the integral 
with respect to 0.^-^ of that product extends into a set of finite measure and 
then is a bounded C°°-function of ui with compact support. Finally, similar 
arguments can be applied to its derivatives. 

The same reasoning can be applied when, instead of I?(]R'^), we consider 
the Schwartz space 5(]R''), whenever the Jacobian Jij) and its derivatives are of 
polynomial growth. 

Since <f>(wi) G X>(M), the integral 

is improper only at ui = when a < < 6. Then, for the study of the 
convergence of this integral it is convenient to obtain an asymptotic development 
of for small values of ui. To this end, let us consider the functional 

depending on the complex parameter A 

{Wl,(j))^[ W^ik)(j){k)dk. (38) 

Jw>Q 

If the C°°-function W is such that the equation W{k) = defines locally a d — 1 
dimensional surface of regular points, that is, for each point fco in the surface 
there exist a neighbourhood V of fco in and a local system of coordinates 
Ml, . . . , ltd such that W{k) = ui for every k &V (for example, when VVF(fco) 7^ 
0), and the same is valid for the equation W{k) = c with c > 0, then the 
distribution is meromorphic with singularities the sequence of simple poles 
H sect. 3.4.2] 

A = -1,-2, ... , -n, . . . 

The residue of the function H38|) at each of these poles can be expressed by 
means of the test function $(ui) defined in formula H37|l . being the residue at 
X = —n equal to 

Res((W^^,0),A = -n) = ^ ^. 

(n — 1)! 

In terms of the distributions 6^'''>{W), we have then 

Res((I^^, 0), A = -n) = tllll (w^), 0) 
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and we can say that the residue of W'^ at the simple pole A = — n is 

Res(M^^, A = -n) = <5("-i) (W) , 

(n — 1 j! 

in complete analogy with the unidimensional case for the distribution .t^. 

Since the behaviour of ^(ui) for ui > e > has not influence over the 
singularities of the integral 



/"OO 

Jo 



the knowledge of these singularities permits us write an asymptotic development 
of for small values of ui. Indeed, in our case [21 sect. 3. 4. 5] 

~ J2 ^^-n {^^''HW),(t)) for ui small. (39) 

From these results it is easy to derive the following 

Proposition 5 Let W{k) = w(fc) — to be a C°° -function, except perhaps for 
some closed set E of singular points with d- dimensional Lehesgue measure zero, 
such that the equation W{k) = defines a d— \ dimensional surface Sq of regular 
points, that is, for each point ka in the surface there exist a neighbourhood V of 
ko in R'' and a local system of coordinates ui,. . . ,Ud such that W{k) = ui for 
every k Cz V (for example, when VW^(fco) 7^ Oj; ond the same is valid for the 
equation W{k) = c with c € (a, b), where M.'^\E = Uf.e(a.b)'^c- Then, given a test 
function (p G I?(R'^), we have: 

(i) When a < < b, the integral 

(— i — ,m)=f -^^dk^ f -^{u,)du, 

is finite if and only if {S{uj{k) — lu), (j){k)) — 0. 
(a) When a < < b, the Cauchy principal value 

(P.P.— ^ ,4>{k)) =Hm/ -4P—dk 



= lim 

e->0 



I — <^(ui)dui+ I — ^{ui)dui 

Ja Ui Ui 



is always finite. 

(Recall that, when ^ [a,b], the integrals in (i) and (ii) are always finite.) 
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Proof: Since G f (R), the asymptotic development of is 

vahd in a neighbourhood of ui = and also we can assume that a and b are 
finite. 

(i) . The unidimensional integral diti converges if and only if for 
any a > we have $(Mi,a;) ~ u" as ui 0. By formula (|39|l . this condition is 
satisfied if and only if {d{uj{k) — lj), (j){k)) = 0. 

(ii) . It is well known ^ theor.1.35] that for a function of the form ^{ui)/ui, 
$ being continuous in a neighbourhood of ui = 0, the integral in the P.P. sense 
exists. □ 



Corollary 6 Proposition\^is satisfied also for every (j) G 5(R'*) if in addition, 
the Jacobian Jip and its derivatives are of polynomial growth. 

Proof: Recall that, if G then also $5(]R'') when the Jacobian J-0 

and its derivatives are of polynomial growth. □ 



Example 7 For the radiative dispersion w(/c) = |fc| — ' take 

ui = uj{k) ~uj,U2^0i, ...,Ud = 9d-i, (40) 
where 9i, . . . , 9d~i are the usual angles in spherical coordinates, we obtain 

being das^_^^^ the Euclidean element of surface for the sphere Sui+u> with centre 
the origin and radius ui + ui. In this case we have 



m 



U){k) — LU 



0('0(w)) dcrs„i- 



dui 



(41) 



Let us put <l>(ui, w) = ^ das^^j^^ 



as 



For w < the integrals of Eq. (|41|l converge because $(ui,a;) is in 
function of mi. 

If w > 0, we are under the hypothesis of proposition[Sland then the integrals 
in Eq. H41(l converge if and only if ((5(w(fc) — uj),(j){k)) = 0, and they converge 
in the P.P. sense for all G V{R'^). 

If w = 0, we cannot apply proposition since the equation W{k) = Lu{k) = 
does not define a regular surface else a singular point fc = 0, but in this case the 
integrals of Eq. converge because $(mi, 0) as mi ^ 0. 

By corollary El the same results are satisfied for cf) G S{R'^). 

Example 8 For Lj{k) = k^ = X]j=i ^j^ applying again the change of variables 
(|40|) . now we have 

^1/1 



1/2 



^'^s r—- = ^ . , das r- 



1+" 
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Then, in this case we obtain 



m 



■ dk 



1 



2-\/"i + ^ 



das 



dui . 



(42) 



Here, <^{ui,uj) — Jg 



da.Q 



and we can apply the same ar- 



guments than in example [7\ to determine the convergence of the integrals in 
(021) • Thus, for u; < the integrals of g^J converge for every e T>{R'^) 
and, on the other hand, when w > these integrals converge if and only if 
{S{u}{k) — uj), 4>{k)) — and they converge in the P.P. sense for all e 2?(M'^). 
By corollary El the same results are satisfied for </« e S{R'^). 



These results are of applicability in the following physical example. 



7 The Hydrogen Atom in the EM Field 

It is well known 3 that in a central potential, caused here by the hydrogen 
nucleus, the bound states of a spinless electron are determined by three quantum 
numbers n, I and m. 

The total or energy quantum number n, whose range of values is n = 
1,2,3,..., +CX), determines the energy E„ of the electron 

1 Zme^ Ze^ 
""~2^^ = "2^' ^ ' 

where m and e are the mass and charge of the electron, oq is the Bohr radius 

ao = ^ = 0.53 X 10"^ cm (44) 

and, for the hydrogen atom, Z = 1. 

The orbital and magnetic quantum numbers, I and m, determine the an- 
gular momentum and the angular momentum along the axis of quantization, 
respectively, and their ranges of values are I — 0, 1, 2, . . . , 7i — 1 and m = 
-1,-1 + 1,. ..,+1. 

In the following q shall denote the position of the electron in the 3-dimensional 
space and we shall assume that the nucleus of the hydrogen atom is fixed at the 
origin. 

In spherical coordinates, the associated total eigenfunctions V'nim are given 

by 

^nU\qlO,4>) ^ Rnl{\q\)Yi.„,{0,(j,) , (45) 

being the radial eigenfunction corresponding to the quantum numbers n and Z, 
Rni, equal to 



RnlM) 



o n 1/2 

2 y {n-l-l)r 

naoj 2n[{n + iy.]^ 



e-l^l/2|5|'L^/+i(|g|), (46) 
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where L'^^^ is the associated Laguerre polynomial 



li!:h\i\)= e 



^ (n- I -1- 3)1(21 + 1 + sV.sl 

s—0 ^ ^ ^ 



(47) 



and the spherical harmonics of order Y;™, are given by 

"2; + 1(;-to)!"1^/^ 



47r {l + my. 



P;"(cos(6')) e*"* 



-1) 



l-\-7n 



1 



d{cos{e)) 



{21 + 1)1 (/-to)! 
4tt (l + m)\ 

sin2'(6') e™"^ 



1/2 



sinl"l(6') (48) 



where P™ denotes the associated Legendre function of the first kind 



(49) 



In the following we shall consider only states of the electron with orbital 
number I = 0. RecaU that yoo(6',0) = (47r)-i/2. 

For Quantum Electrodynamics, the expansion of the EM vector potential in 
Fourier integral is 

A{q) - [ dkg{k)e'''-iak, 



where, as until now, q denotes the position of the electron in the 3-dimcnsional 
space, and fc G M'^ corresponds to momentum coordinates. From this expansion 
we obtain the response terms 



D{k) 



|fc|l/2^ 



(50) 



being the cutoff or form factor g{k) — \k\ 



-1/2 



Let us introduce some notation. For the positive Bohr frequencies we shall 
write, see Eq. (p^ . 



uj ■— Fj — E 



(1. _ 1. 

2ao V TO-^ 



, m > n , 



and for the matrix elements of the operators D{k) we shall put 

(k) (V'm00,-D(fc)V'n00) 

- {Rrno{\q\)Yoo{0, 4'),D{k) RM\q\)yoo{0, </>)) ■ 



(51) 



(52) 



17 



Now suppose that for the positive Bohr frequency tOmn there exists a unique 
pair of energy levels Sm — E„i, £n = En in Spec Hs such that oJmn = Em ~ En-^ 
Then, the operator Di^^^^ defined by Eq. (|^ can be written as 

(k) |'0nOo)(V'mOo| 

and the imaginary part of the corresponding generalized susceptivity factor 'y^™" 
is of the form, see Eq. 

/m(7^-) = -P.P. / kl^^IiSrffc. (53) 

In what follows we shall study the generalized Hilbert transform of the right 
hand side of Eq. (|53|l in terms of the cutoff function g. We restrict our attention 
to the dispersion function uj{k) = \k\ of exampled For u!{k) = fc^ the results 
are very similar. 

7.1 Cutoff Functions of the form (z/ > 0). 

\kr 



When the response terms under consideration are of the form 

Dik)^^^, iy>0, (54) 



the matrix element corresponding to the states of the electron with energy 
numbers m and n and orbital number ^ = is 

9nUk) = jjj^ E ((TTw " i^-myj ' ^^^^ 



where 



^ T\s A min{ri-l,s-2} , -, \ / i \ / 

(-1) 4 / n~l \ I m-1 \ / s 



u ^ ' Q=max{0,s-m-l} ^ ' ^ ' ^ 

(56) 

The change of variables (|4()|l gives us a more convenient expression of the 
right hand side of Eq. (|53|) . For the dispersion function ti;(fc) — \k\ we have 

P.P. / lg""WI' rffc = p,p 



m+n 

mn 



s=2 



w(fc)-wm„ ■ ■ ui (ui + Wm„)2'' 



(1 + 2(ui + W„„))' (1 - i(ui + iOm-a))'' 



We can already determine the values of the parameter v for which the imag- 
inary part of the generalized susceptivity factor ■y"'"" exist. 

^This is not the case for every positive Bohr frequency LOmn of the hydrogen atom, but 
by means of a httle perturbation the hydrogen atom becomes a system verifying this generic 
assumption. 
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Proposition 9 Let us consider the interaction of a spinless electron in the 
hydrogen atom with an electromagnetic field, in which the response terms are of 
the form 

^(fc)-^, ^>0. (58) 

Then, for the dispersion functions uj{k) — \k\ (and also for Loik) — kF'), the 
Cauchy Principal Values 

pp f lg""WI' dk, m>n EN, (59) 
are finite if and only if v < 3/2. 

7.2 Cutoff Functions of the form g{\k\). 

Now, let us consider response terms of the form 

i?(A:)=g(|fc|)e*«, (60) 

being the cutoff function g, at first, a medible fmiction of \k\. 

In this case the matrix element corresponding to the states of the electron 
with energy numbers m and n and orbital number ^ = is 

*5(|fc|) v-^^m„ ( 1 1 



where the coefficient C™" is given by Eq. H56|l . 

For the dispersion function u){k) = the change of variables (|40|l gives us 
now 

pp I dfc-p.p. r ^\9{ni+^n.n)\\ ,^ 

m+n / \ 2 (62) 

\ ^riin I 



From the last expression we can deduce sufficient conditions on the form 
factor g in order that the imaginary part of the generalized susceptivity factor 
7""" exist. 

Proposition 10 Let us consider the interaction of a spinless electron in the 
hydrogen atom with an electromagnetic field, in which the response terms are of 
the form 

i?(fc)=g(|fc|)e*^ (63) 
Then, for the dispersion function uj{k) = \k\, the Cauchy Principal Values 

P.P. / lw(fc)P m>neN, (64) 

Jr3 U)(k) - LOmn 

are finite if the function g verify the following conditions: 
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(al) \g{ui +ijJ„i„)p/uf G L^(\b,(X))) , with respect to ui, for some b > 0; for 
example, if \g{ui + u!mn)\ — u\ as ui go, for v < 3/2; 

(a2) the integrand of the right hand side of Eq. \6S\) . in a neighbourhood of 
Ml — 0, is the sum of an antisymmetric function f\ and a symmetric 
function f2 such that fii^i) dui is finite for some e > 0; for example, 
when g is bounded in a neighbourhood ofujmn (or, equivalently, g{ui+U!mn) 
is bounded for ui in a neighbourhood of ); 

(aS) \g{ui+LUrnn)\'^{ui+U!mn)^ G ( [- , - + e]) with rcspcct toui, for 
some e < uj^n; for example, if\g{ui+uimn)\ — {ui+iOmnY a* "i ^ -^mn, 
for V > —3; 

(a4) g{ui + ijJmn) S L'^iy^t^mn + ^e] U [c, &]), with respect to ui, for some 
< e < ujmn and any finite b > e. 

Note that condition (a2) implies some special behaviour of the cutoff function 
g in a neighbourhood of the resonance surface ui — uj{k) — cOmn — 0. 
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